Where are the Fermi Lines Coming From? 



Kanishka Rao 1 and Daniel Whiteson 1 
1 Dept. of Physics & Astronomy, UC Irvine, Irvine, CA 

We estimate the spatial locations of sources of the the observed features in the Fermi-LAT photon 
spectrum at E-, = 110 and _E 7 = 130 GeV. We determine whether they are consistent with emission 
from a single source, as would be expected in their interpretation as 77 and 7Z lines from dark 
matter annhiliation, as well as whether they are consistent with a dark matter halo positioned 
at the center of the galaxy. We take advantage of the per-photon measured incident angle in 
reconstructing the line features. In addition, we use a data-driven background model rather than 
making the assumption of a feature-less background. We localize the sources of the features at 110 
and 130 GeV. Assuming an Einasto (NFW) density model we find the 130 GeV line to be offset from 
the galactic center by 285 (280) pc, the 110 GeV line by 60 (30) pc with a large relative separation of 
220 (240) pc. However, we find this displacement of each source from the galactic center, as well as 
their relative displacement to be statistically consistent with a single Einasto or NFW dark matter 
halo at the center of the galaxy. 
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INTRODUCTION 

The search for the particle nature of dark matter is 
one of the oustanding open questions of modern physics. 
A broad program of research with a variety of promis- 
ing and complementary approaches has become a major 
piece of experimental particle physics research. This in- 
cludes searches for direct production at colliders, searches 
for dark matter interactions with standard model parti- 
cles in large quiet underground detectors, and searches 
for signals of dark matter annhiliation into standard 
model particles in regions of the galaxy with large dark 
matter density. 

A clear signal of dark matter annihilation may be car- 
ried by gamma rays traveling to Earth from regions in 
the galaxy with high dark-matter density. As the pho- 
tons do not typically scatter after their production, the 
photon energy and direction are powerful handles for un- 
derstanding the mechanism of dark matter annihiliation 
into standard model particles. 

One mechanism is annihilation resulting in quarks, 
which would hadronize and yield ir° particles which in 
turn produce photons. The broad spectrum of such a pro- 
cess would give fairly low energy photons {E 1 < 50 GeV) 
and may be difficult to distinguish from other sources. 

A clearer signature may appear from annihilation di- 
rectly into two-body final states including a photon, 
though the rate would be smaller than continuum emis- 
sion due to loop supression. Rather than yielding a broad 
energy spectrum, this process would produce a photon 
with a well-defined energy (a "line" ). This makes a search 
for lines in the photon spectrum an important component 
of the dark matter program using Fermi-LAT data [IH2] . 

Recently, observation of a feature with high local sta- 
tistical significance at £L = 130 GeV was reported [H[S], 



with a source location 1-2 degrees away from the galactic 
center jB]. Follow-up analyses suggested the possibility 
of two features, consistent with the spacing of lines ex- 
pected from 77 and 7Z processes [7J |S] 

In this paper, we study in detail the question of the 
location of the source of the photons in the features. 
We confirm the location of the photons in the feature at 
E 7 = 130 GeV, and for the first time, locate the source 
of the photons in the feature at E 1 = 110 GeV. We study 
whether the two are - individually or collectively - con- 
sistent with emission from a single source, as would be 
expected if they represent the 77 and 7Z processes, as 
well as whether that source is consistent with a dark mat- 
ter halo at the galactic center. In addition, we introduce 
a new approach to analyzing the data which uses the 
measured angles of the individual photons, rather than 
the median photon angles. 

ENERGY SPECTRUM ANALYSIS 

A source of photons with a very narrow range of en- 
ergies would appear in the spectrum of observed Fermi- 
LAT photons as a peak, due to finite energy resolution. 
Previous analyses [H [7J [5] have used the Fermi-LAT en- 
ergy dispersion formula [3] to build a probability density 
function (pdf) for the reconstructed energies of photons 
from a line, and fit it to the observed spectrum. 

Given the small size of the dataset, we seek to use 
the maximum available information about each photon 
to extract insight into the observed features. The energy 
resolution, for example, has a significant dependence on 
two other photon characteristics: 6, the incident angle 
relative to a line normal to the LAT face, and the photon 
'type' which indicates whether the conversion happened 
in the front or back layers of the tracker [TIT] ■ Figure [I] 
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shows the reconstructed energy pdf for various choices of 
9 and type. This per-photon information is very perti- 
nant to the analysis of this spectrum. If the photons in 
the peak were all well-measured, it would enhance the sig- 
nificance of the peak; conversely, if they were all poorly- 
measured, it would degrade the significance. 

Our approach is to reconstruct the distribution of 
i?* ruc , the true energy of photons striking the LAT, rather 
than the distribution of E^ eco , the reconstructed energy 
of each photon. The former distribution can account for 
the per-photon resolution, while the latter treats each 
photon equally. The reconstructed energy is related to 
the true energy by the pdf f (E Ieco \E tIue , 9, type). We 
use Bayes theorem to invert this relation and calculate 
/(£ truo |£: rcco ,(9,type) where P(E Icco ^ is taken from the 
full-sky photon energy spectrum. Other choices of pri- 
ors have negliglible effect, as the energy dispersion pdf is 
sharply peaked, see Fig[T] 

The distribution of £"' ruc is then built as a kernel 
density estimate [TQ with f (E true \E reco , 9, type) as the 
choice of kernel. Other approaches |6] have used 
adaptive Gaussian kernels, but our physically motivated 
choice of kernel reflects the asymmetric shape of the pdf 
and uses - for the first time - the per-photon measured in- 
cident angle rather than the median angle of the dataset. 

We use the publicly available Fermi-LAT data collected 
over 3.7 years. We use ultraclean photons with zenith 
angle < 100 to veto photons from Earth as well as stan- 
dard quality requirements [12]. We select a three-degree 
region around the galactic center. 

The observed photon counts in the three-degree region 
are shown in Fig.pl The distribution of i?' rue is shown in 
Fig. [3j with several choices of angle-dependence used to 
demonstrate the impact of this information. Two choices 
demonstrates the range of possibilities, if all photons had 
the maximum or minimum possible energy resolutions. 
This demonstrates the potential importance of this in- 
formation. The spectrum which uses the measured per- 
photon angles has peaks at 110 and 130 GeV which are 
slightly degraded relative to a spectrum which uses the 
median photon angle. As the distribution of photons in 
the peak region and the background region have similar 
distributions of cos(#) and conversion type [13], there is 
not currently a large impact from using the per-photon 
information, but it remains a useful method to visualize 
the spectrum. 
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FIG. 1: Probability density functions which describe the 
energy dispersion of the Fermi-LAT. Top, the probability 
density for reconstructed photon energies for a photon with 
E* IUC = 130 GeV. Bottom, the probability density for true 
photon energies for a photon with E™ c ° = 130 GeV. 
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FIG. 2: Reconstructed photon energies Elf c ° from Fermi- 
LAT in a three-degree region surrounding the galactic center. 



LOCATING THE FEATURES 
Technique 

We search for the most likely source position using a 
maximal likelihood fit of the observed photons to the ex- 
pected spatial and energy distributions from dark matter 
halos. The unbinned likelihood is 



m,b)=ii 



n b g fbg(Ei, 9i) + n sig f sig (Ei, li,bi\9i, I, b) 



n hg 



i-sig 



(1) 

where (I, b) are the galactic coordinates of the cen- 
ter of the dark matter halo, n hg = and the 
background and signal normalizations, the index i runs 
over observed photon energy and spatial parameters 
(Ei,9i,li,bi) within the energy window surrounding the 
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FIG. 3: Relative likelihood of true photon energies, _E* rue , 
built from the reconstructed energies of Fig [2] and the pdf 
/(S truc |£; rcco ,f5,type), see text for details. The impact of 
different assumptions about the photon incident angle 8 is 
shown. 



feature(s) of interest, and f s i g (Ei, Zj, bi\6i, I, b) is the dark 
matter halo pdf as a function of photon energy Ei, in- 
cidence angle 0j at spatial position (li 7 bi) given a dark 
matter halo centered at (l,b). The pdf / s ; g accounts for 
per-photon spatial [14| and energy [S] resolution in the 
same spirit as the Bayesian unfolding described above; 
here we do not need to explicitly unfold, as we have in- 
cluded the per-photon information in the unbinncd like- 
lihood. 

The background pdf /b g is built from a data-driven 
model, with the energy dependence coming from photons 
outside the three-degree region surrounding the galac- 
tic center, see Figure [3] We find this to be reason- 
ably consistent with a power law model used in previ- 
ous analyses @HE]. We use energy windows of E^' uc = 
[105, 115], [125, 135] for the location of the feature at 110 
GeV and 130 GeV, respectively, or both windows for the 
combined features. In the case of the two-feature analy- 
sis (Ej = 110, 130 GeV), there are two / s j g terms in the 
likelihood, one for each peak. We use normalizations of 
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6 and 



14 for the £L 



110 and 130 GeV 



features, respectively. 

The dark matter halo pdfs / s ; g are derived from either 
NFW [T5] or Einasto [TB] halo profiles; the pdf in (7, b) is 
calculated via the linc-intcgral of the square of the dark 
matter density [17 . We use Q£ = 0.17 for the Einasto 
model and r s = 20 kpc for both models. 



Results 

The position of the most likely values of the halo cen- 
ters are given in Table [I] for both DM profiles and each of 
the features as well as the combined spectrum. In both 
cases, the feature at 130 GeV appears to be displaced 
from the galactic center, as previously reported, but the 



FIG. 4: The data-driven background model (points) con- 
structed from photons outside the three-degree region sur- 
rounding the galactic center, compared to a power law with 
exponent -2.6. 



110 GeV feature appears to be centered at the galactic 
origin. 
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FIG. 5: Most likely positions in galactic coordinates (l,b) of 
an NFW (top) or Einasto (botton) dark matter halo, max- 
imized seperately for the features at 130 GeV, 110 GeV, or 
the combined features. Also shown are la and 2a uncertainty 
regions, see text for details. 

Figure [5] shows the most likely positions for each of 
the dark matter halo profiles and each of the energy 
spectrum features, as well as uncertainty regions. The 
uncertainty regions are calculated in a frequentist man- 
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TABLE I: Most likely positions in galactic coordinates (l,b) 
of an NFW or Einasto dark matter halo, maximized seper- 
ately for the features at 130 GeV, 110 GeV, or the combined 
features. 



130 GeV 110 GeV Combined 
"NFW (-1.5,-1.1) (-0.2,0.1) (-0.4,-0.2) 
Einasto (-1.5,-1.2) (-0.4,-0.2) (-0.2,0.1) 



ner, using simulated experiments to determine the AL 
threshould which would contain the true (I, b) position in 
68% (for la) or 95% (for 2a) of the cases. 

If the two features were both due to dark matter 
annhililation processes, one would expect them to be co- 
located. Figure [6] shows the unfolded E^ uc spectrum in 
a 1-degree circle surrounding each of the features; in the 
spectrum which highlights the 130 GeV feature, the 110 
GeV feature is suppressed, and vice versa. Similarly, Fig- 
ure [7] shows the unfolded £'* ruo spectrum in the la region 
surrounding each feature. 



Best NFW 1° region 
Best Einasto 1 ° regio 



.+++4 



Bes! NFW 1° region 
Besl Einasto 1° region 



FIG. 6: Spectrum in a region within 1 degree of the maxi- 
mum likelihood for 130 (left) and 110 (right) features assum- 
ing an Einasto or NFW source. 
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FIG. 7: Spectrum in a region within la of the maximum 
likelihood for 130 (left) and 110 (right) features assuming an 
Einasto or NFW source. 

This suggests that the two features of the energy 
spectrum may have spatially distinct sources, and these 
sources may themselves be offset from the galactic cen- 
ter. To evaluate the significance of the sources' displace- 
ment from the galacter center and from each other a more 
rigorous analysis which accounts for uncertainties is per- 
formed below. 



Interpretation 

If the observed features are both due to dark matter 
annhilation, then their locations in the sky depend only 
on the dark matter density and not the nature of the final 
state, 77 or 7Z, and so should yield locations consistent 
with a single source. As shown above, we find that the 
sources of the E 1 = 110 and 130 GeV features are seper- 
ated appreciably from each other, and that the E~ = 130 
feature is displaced from the galactic center. 

Definitive conclusions, however, require analysis of the 
uncertainties involved, including the expected fluctua- 
tions of the measured source locations given the small 
number of photons in each feature. We seek to answer 
two questions: 

• Are the two features consistent with emission from 
a single source? 

• Are the features consistent with emission from a 
dark matter halo at the galactic center? 

We answer these questions by comparing hypothe- 
ses using likelihood ratios. For example, we can probe 
whether the feature at E 1 = 130 GeV is consistent with 
emission from the galactic center by evaluating 



-2 log 



L(l = l,b = b) 
L{1 = 0,6 = 0) 



(2) 



where (I, b) is the location which maximizes L and (0, 0) 
serves as the null hypothesis. When (I, b) is close to (0, 0), 
the log of the likelihood ratio q approaches zero; when it 
is far from the origin, it becomes negative. We compare 
the measured value of q to the expected distribution in 
simulated experiments from an NFW or Einasto profile 
with a line at E 1 — 130 GeV placed at (0,0), see Fig- 
ures^ a) and [9^ a). While the value of q indicates that 
the maximal position is displaced from the origin, it is 
consistent with the expected distribution of q values when 
the DM halo is at the origin; the p- values = 0.18 and 0.14 
for Einasto or NFW, respectively. Distributions of q for 
DM halos offset from the origin are also shown. 

We perform this analysis for the feature at E 7 = 110 
GeV as well, see Figures (8^b) and (9^b). As this feature 
appears to be located close to (0,0), the q values are 
close to zero and entirely consistent with a DM halo at 
the galactic center. 

We repeat the analysis for the combined features at 
110 and 130 GeV, where we constrain them to have the 
same location, see Figures [|Jc) and |9jc). Again, the 
q values are consistent with a DM halo at the galactic 
center. 

To address the question of whether the two features are 
more likely to be from a single source (consistent with a 
DM interpretation) or from two spatially distinct sources, 
we construct a likelihood ratio 
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-2 log 



Li3o(l,i>) x L 110 (l,b) 

il3o(^130> ^130 ) x illo(^110j ^lio) 



(3) 



where in one case we find the single location (l,b) which 
maximizes L for both features, and in the other we allow 
the two features to find individual maxima: (lno 5 ^no) 
and (/130, &i3o)- This q will be zero if the individually op- 
timized (and possibly separated) locations are consistent 
with the most likely single joint location, and positive 
if individually optimized (and possibly separated) loca- 
tions give a better description of the data. Again we 
compare to the expected distribution of the likelihood 
ratio in simulated experiments in which the two sources 
are co-located and find that the observed q is consistent 
with a co-located source, see Figures [§d) and gd). 
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FIG. 9: The log of likelihood ratios comparing the likelihood 
of an NFW DM halo placed at the most likely position (I, b) or 
at the galactic center (I = 0, b — 0), in the case of the 130 GeV 
feature (a), the 110 GeV feature (b), or the combined features 
(c). Also shown is the likelihood ratio between two co- located 
sources or sources which are seperated (d). In each case, we 
show the expected distribution in simulated experiments. 
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FIG. 8: The log of likelihood ratios comparing the likelihood 
of a Einasto DM halo placed at the most likely position (I, b) 
or at the galactic center (I = 0, 6 = 0), in the case of the 130 
GeV feature (a), the 110 GeV feature (b), or the combined 
features (c). Also shown is the likelihood ratio between two 
co-located sources or sources which are seperated (d). In each 
case, we show the expected distribution in simulated experi- 
ments. 



different from the standard approach of assuming the me- 
dian angle and therefore resolution. 

We have presented an analysis of the spatial locations 
of the sources of the observed features in the Fermi-LAT 
energy spectrum at E 1 = 110 and 130 GeV. We find: 

• the 130 GeV photons come from a location dis- 
placed from the galactic center, but consistent with 
an Einasto or NFW dark matter halo at the center. 
This is consistent with previous work [TT] . 

• the 110 GeV photons come from a location close to 
the galactic center, and consistent with an Einasto 
or NFW dark matter halo at the center. 

• the two features at 110 and 130 GeV photons are 
consistent with a single DM halo at the galactic 
center, with either an Einasto or NFW density pro- 
file. 



CONCLUSION 

We haved demonstrated via Bayesian unfolding the po- 
tential importance of including the per-photon energy 
resolution, though for this dataset the results are little 



Some have suggested that the displacement of the 
source of the E 1 — 130 GeV feature indicates that DM 
halos may be displaced from the galactic center [T5] , but 
we find that the current data are consistent with a DM 
halo from the center. 
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